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G\ ■ Abstract 

i-^ \ Depending on their masses, Supersymmetric particles can affect various mea- 

[t \ surements in Z decay. Among these are the total width (or consequent ex- 

tracted value of a s ), enhancement or suppression of various flavors, and left- 



m 



right and forward-backward asymmetries. The latter depend on squark mass 
splittings and are, therefore, a possible test of the Supergravity related predic- 



cn 

>. 

t — ■ , tions. We calculate leading order corrections for these quantities considering 

in particular the case of light photino and gluino where the SUSY effects are 

^\j ' enhanced. In this limit the effect on a s is appreciable, the effect on R^ is 

small, and the effect on the asymmetries is extremely small. 

PACS numbers: 12.60.Jv, 12.38.Bx, 13.35.Dx, 13.38.Dg 

^ ■ 

a,: 
a: 

•i-H 

X 



Typeset using REVTgX 



*lclavell@ualvm. ua.edu 

fpcoulter@ualvm.ua.edu 

flevan@gluino.ph.ua.edu 



1. INTRODUCTION 

In recent years several slight anomalies in Z decay have encouraged attention to the 
possibility of supersymmetric particles in the Z region and below. These are 1) an apparent 
value of the strong coupling at the Z greater than expected from the extrapolation of low 
energy measurements, 2) an enhanced value of the Z branching ratio into b quarks relative 
to the Standard Model (SM) expectations, and 3) deviations from the Standard Model 
predictions for the left-right asymmetry in polarized electron-positron annihilation and for 
various forward backward asymmetries. Even if these discrepancies shrink under further 
analysis as may now be the case with the Rb anomaly, it is useful to consider what constraints 
they impose on SUSY masses. In this note we treat the strong supersymmetric (SUSY) 
effects on these quantities from the graph of Fig. 1 involving a virtual squark-squark-gluino 
triangle correction to the Z coupling to quark-antiquark final states. We assume that, while 
the gluino might be light relative to the Z mass, the squarks are all above the Z so there 
are no effects from production of on-shell squarks. The graph of Fig. 1 has been considered 
earlier fl|-§J. These works for the most part ignore sfermion mass splittings crucial to some 
of the anomalies discussed above and, in any case, do not consider explicit mass splitting 
schemes currently under discussion. The present work, extending the results of these authors 
to include these effects, is divided as follows. In section 2 we briefly discuss the experimental 
situation relative to the above-mentioned anomalies. In section 3 we present the calculation 
of the graph in Fig. 1, allowing for squark mass splittings. Results and conclusions are given 
in section 4. 

2. ANOMALIES IN Z DECAY 

Since 1992, various authors have noted that the apparent value of the strong coupling 
constant at the Z is greater than the value expected in the Standard Model from extrapo- 
lation from lower energies. This was in fact one of the initial proposed positive indications 
0-0] of a light gluino since such a particle would slow the running of a s . The effect has been 
clouded by disagreement over the value of this coupling at both the low and high energy 
regions but, clearly, the greater is the value of a s (M z ) and the lower is its value at lower 
energies, the more likely one is to be interested in the light gluino possibility. The fitted 
values quoted from Z decay measurements have ranged from a low of 0.108 ± .005 || to a 
high of 0.141 ± .017 ||. Current values are typified [|K| by that from the Z hadronic width. 

a s (M z ) apparcnt = 0.123 ± 0.005 (2.1) 

Here we have labeled the fitted value as " apparent" since the value quoted is from a Standard 
Model analysis and there could be extra (e.g., SUSY) contributions such as that of Fig. 1. We 
assume that the actual value of a s (Mz) is that extrapolated from low energy. Although these 
numbers are somewhat dependent on the analysis performed, they tend to be systematically 



below the value in Eq. ( |2.1| ). For definiteness we will note the values from [|I(J in the SM or 
heavy gluino (HG) case as obtained from an T decay analysis 



a s (M z ) bM = 0.108 ±.010 (2.2) 




FIG. 1. Leading order virtual squark-squark-gluino correction to hadronic Z width. Solid and 
dashed lines denote quarks and squarks correspondingly. Wavy line denotes gluino 



and from [PI] in the light gluino (LG) case 

a s (M z ) LG = 0.114 ±.005. 
Thus the difference between the actual a s (M z ) and the apparent value is 



5a. 



5a s = -0.009 ± .007 (LG case) 
-0.015 ± .011 (SM or HG case) 



(2.3) 



(2.4) 
(2.5) 



The SM value of the coupling constant adopted here is in agreement with that found from 
QCD sum rules |I^1 , although these authors quote large errors. Other particular analyses, 
notably that from r decay, could significantly reduce the discrepancy (or even change it's 
sign) but we feel it fair to say that a large body of low energy data prefers negative 5a s . 



The t analysis in particular has been criticized |13j as being especially vulnerable to non- 
perturbative effects. In the LG case there is a small positive contribution to 5a s from gluino 
contributions alone to Z decay [14]] amounting to 5a s ~ +.002. Related perhaps to the 
a s anomaly, is a possible excess of b quark production in Z decay which has stimulated 
significant interest in possible manifestations of SUSY below the Z. In particular it has been 
noted that an appreciable enhancement of b production could be a signal for charginos and 
stop quarks below the Z PJI5|. Such a non-standard-model contribution to Z decay might 
also explain the apparent enhancement in a s at the Z relative to the value expected from 
extrapolation from lower energies. If the only non-SM contribution to the Z width was in 



the bb channel it would be related to 5a s by [|16[ 



5 a „ = irSR, 



1 + a s /ir 
l-Rb '' 



where 



Ri> 



T{Z -* bb) 
Y{Z —> hadrons) 



(2.6) 



(2.7) 



and 5Rb is the deviation of this quantity from its standard model value. However, due to 
new LEP results above the Z, the likelihood of low mass charginos (at least in the standard 



SUSY picture) has diminished. In addition the b anomaly has largely disappeared under 



reanalysis. The latest experimental analysis yields |17 



5R b = .0009 ±.0011 (2.8) 

This value is consistent with zero and it is interesting to see how severely it constrains 
contributions beyond the Standard Model and, in particular, whether SUSY effects could 
lead to an enhanced apparent a s without enhancing the b decay mode of the Z. With 
sufficiently light chargino and stop quarks the chargino-stop virtual correction to the Z 
vertex dominates the SUSY contribution to the Z partial widths and can lead to a SRb as 
great as .0018. In the light gaugino scenario, however, with mi/2 = A = 0, the stop quarks 
are expected to be above the top [TH] and this contribution to 5Rb drops below .0005. In 



this case, the virtual gluino correction could become the dominant SUSY effect. In this 
article we reconsider the gluino (and photino) corrections with special attention to scheme 
dependence and dependence on the sfermion mass splitting which have not been part of 
prior investigations. 

We also consider the left-right and forward-backward asymmetries. These are related 
to the effective coupling of the Z to vector and axial vector currents of particular fermion 
species. One defines 

A(f) 2g v (f)g A (f) 

The left-right asymmetry is defined in terms of the total annihilation cross sections for 
left-handed and right-handed polarized electrons. 

A LR = °J^-^ = A(e). (2.10) 

Vl + ctr 

The Forward Backward asymmetry for a particular flavor f is 

a m _ ^(cosg>0)-g(cosfl<0) 
FB{J) a(cos6>0) + a(cos6<0)' { ' ' 

For an unpolarized beam this is given by 

A FB (f) = \A{e)A{f). (2.12) 

Current asymmetry measurements can be summarized by experimental determinations of 
the difference between the A(f) and their standard model predictions fT0|| . 



6A{e) = 0.011 ± .005 

8A{b) = -.093 ± .053 (2.13) 

SA(c) = -.061 ± .090 

Even if, as we find, the SUSY contributions cannot account for discrepancies at this level, 
it is useful to record here the SUSY prediction for comparison with future more accurate 
measurements. 



3. SUSY QCD VERTEX CORRECTION 

The Born term matrix element for Z decay into a fermion-antifermion pair in the Stan- 
dard Model takes the form 

m° 7 = -A^^iphM - ghsMq - p), (3.1) 

J J Sill Au 

where 

g v = T 3 -2Qsin 2 e , g° A = T 3 . (3.2) 

In terms of these parameters the zeroth order Z width into a fermion anti-fermion pair, 
neglecting the fermion mass, is 

T )j = -^Y 6 M M) 2 + (9a)\ (3-3) 

A second quantity of some interest is 

since the forward-backward asymmetry of fermion species / is given by eq. ( [2.12|) , and 
the left-right asymmetry of the total hadronic cross section with a longitudinally polarized 
electron beam is Apr = A°(e). The existing discrepancies between these measured quantities 
and the Standard Model predictions make a study of the supersymmetric contributions of 
interest. The lowest order SUSY QCD corrections to the Z-fermion vertex arise from the 
squark-squark-gluino triangle graph of Fig. 1. The corresponding matrix element is 



5M z^n = —^ e Mph»( 5 9v - SghsMv - p). ( 3 -5) 

J J Sill Zu 



where 



Sgy = [g v (Il + Ir) + 9WI - IrWf^, (3.6) 

S9A = [g A (i r L + i R ) + g°v(i r L - Ir)}Cf^, (3.7) 

and Cp is the eigenvalue of the quadratic Casimir operator of the gauge group in the fun- 
damental representation (for SU(3), Cp = 4/3). The I r L R are the contributions from left 
and right squarks corresponding to the graph of Fig. 1. The r indicates that the ultravio- 
let renormalization has been made. Adding the SUSY correction to the Born term matrix 
element is equivalent to replacing gy A by 

9v,A = gy,A + $9v,a- (3-8) 

These lead to the 0(a s ) equations 



gUf) + g 2 Af) = U(f)) 2 + (9°Af)) 2 } (i + c F ^ 



(r L + r R ) + A°(f)(ii-r R ) 



(3.9) 



A(f) = A°(f) + C F ^-[ r,-v 



2n 



R 



(3.10) 



Using dimensional regularization [19[ with the dimension of space - time d= 4 — 2e and a 
standard Feynman parametrization, the unrenormalized Il,r can be written as follows. 



or 



where 



2-d,,2e 



l L,R 



16i(27r) 2 -V 
d-A 



I d d k I dxidx 2 0(l — Xi — x 2 ) 



C 



L.R 



{Cl,r - fc 2 ) 3 



I L . R = r( e )(47r)^ 1 (%^)-«r/.r l r/r^i I - ,r,. .v 2 ) 



/'• 



(3.11) 



c 



L,R 



[i x 2 + m LR (l - x 2 ) - Q x 1 (l - x\ - x 2 ). 



(3.12) 



Here \x is the gluino mass and mj,^ are the left and right squark masses of the same flavor 
as the final state fermions. Q 2 is the center-of-mass energy squared (M| at the Z pole) and 
\i is t'Hooft's unit of mass fill . 

Adding SUSY self energy corrections on the external quark legs has the effect of sub- 
tracting from the result of Fig.l its value at Q 2 = @]. Then 



L L,R 



o 



dxidx 2 9(l — X\ — x 2 ) In 



Q 2 xi{\ - xi - x 2 ) 



/j, 2 x 2 + m z LR (l -x 2 ) 



(3.13) 



This function has the series expression 

oo r^ln 



jCWZ 
l L,R 



E 

n=l 



Q" 



n 



-B(n+l,n + l) 



dx 



i-x 2 y> 



2n+l 



2 7T2 



\fi x 2 + m 2 LR {l - x 2 )\ 



(3.14) 



Two simple cases can be studied: 



r ,~ m v-m Q 2 \ B(n+l,n+l) 1 / Q 2 



n=l^ m L,R< 



n(n + 2) 



m 



L,R 



J_(_& 



240 \m 2 



L,R 



+ 



and 



h, R (fi = ™l,r) = E ( ip- ) 



Q 2 \ n B(n+l,n + l) 1 Q 2 



n=i^i,fl- 



n(2n + 2) 



24 m 2 



+ 



(3.15) 



(3.16) 



In each case the series are rapidly converging and the first term gives about 85% of the full 
result even down to fh 2 L R = Q 2 . One can see from these results that although I falls rapidly 
with increasing squark mass it is much less sensitive to the gluino mass. Since an ultra-light 
gluino is still a possibility |2"0 ], in the following section we concentrate on the numerical 
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results for the \x = case. The one-loop SUSY QCD corrected decay rates and asymmetries 
into // are to 0(a s ) 



n 



ff sin 2 26 



M z (g 2 v + g 2 A ) = T° f7 



-C F ^ + C F 

4 71 



Mf) 



29V9A A o {f) + Cp 



9v + 9a 



a s 


(H 


+ Ir 


IT 


{ 


2 


a s 


II- 


jr 
1 R 



A°(f) 



II 



l R 



71 



(3.17) 



(3.18) 



In T r-7 we have included the QCD order a s correction. The total hadronic decay rate to this 
order, obtained by summing the partial rates over /, can be written as follows. 



tot 



r° 



n 



apparent 



7T 



(3.19) 



where the "apparent" a s , the one that would result from a Standard Model analysis, is 



5a, 



Os 



^apparent = ff ^^ 



(3.20) 



where 



s f = c, 



a s r /7 



27cT 



[I r L(f) + I R (f) + A (f)(I r L (f)-I R (f))]. 



(3.21) 



tot 



r° ot being the zeroth order total hadronic decay width. 



The contribution to Rh is 



5R b = 5 b -R° b J25 f . 

f 



(3.22) 



with the analogous expression for 5R C . The Il,r depend on the masses of the SUSY partners 
of the left and right handed quarks of flavor f as well as on the gluino mass. The squark 
masses, in the SUGRA scheme are given by 



z'2 



m 



j L = m 2 + mj + Mz cos 2(3(T 3J - Q f sin" 



mi R = mQ + mj + M z Q f cos 2f3 sin 2 9 



(3.23) 
(3.24) 



where mo is the universal scalar mass. Neglecting the quark masses, and M z , the Il,r 
become independent of flavor. In this case the lowest order contributions to 5R b and to the 
asymmetries vanish although there is still a significant effect on the Z total width. 

4. RESULTS AND CONCLUSIONS 

Since the experimental evidence for discrepancies from the Standard Model are so slight 
(two standard deviations or less), at this point our results should be interpreted as predictions 
of low energy SUSY rather than as suggestions from the data. Our main point is that the 
graph of Fig. 1 can be effective in enhancing the hadronic decay rate of the Z without 
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TABLE I. Results for different schemes. 






/ 


rh L 


fflR 


r° 


5T x 10 4 


SR f x 10 5 


A 


5A x 10 4 


oA oxp 


c 

b 


100 
113 


103 
107 


.287 
.370 


6.73 

6.78 


+ 5.78 
-3.85 


.667 
.935 


0.826 
-1.06 


.06±.09 
.10±.05 



enhancing the relative contribution from b quarks. The variations of SRb and Sa s are 
shown in Figs. 2 and 3 respectively as functions of the universal scalar mass m , using 
a s (Mz) = 0.113, tan/3 = 1.6 and sin 2 # = 0.2317. The contribution to 5Rb and to 5a s 
are negative and slowly varying. The SUSY QCD contribution to 5Rb therefore tends to 
cancel the chargino-stop contribution to 5Rb while still leaving an appreciable negative 
contribution to 5a s consistent with the experimental indications discussed in section 2. 
The contributions to the forward-backward asymmetries behave similarly although they are 
orders of magnitude smaller than the current experimental errors. The SUSY predictions are 
very close to those of the Standard Model and are two standard deviations from the current 
experimental results. In table 1 we summarize the results in the case m = 106GeV and 
jl ~ as suggested in recent phenomenological studies of the Fermilab jet transverse energy 
distributions, scaling ratio, and top quark events PIUI8||. SA exp is defined as A exp — A° and 5T 



is the SUSY correction to the Z width from eq. (|3.17| ). The left-right asymmetry is obtained 
by replacing Cpa s by a in eq. ( |3.18| ). Since the Standard Model result is suppressed by a 
factor of 4 sin 2 9—1 one could expect the experiment to be sensitive to the SUSY contribution. 
However in the SUGRA related model for sfermion mass splittings, the selectron mass 
splitting is also suppressed by the same factor. This fact plus the proportionality to the 
fine structure constant makes this experiment very insensitive to SUSY. 

The analysis presented here can be readily extended to the higher Q 2 values of LEP II. 
The early results from CERN are consistent with the Standard Model but, in view of the 
large fluctuations seen in the early results from LEP I discussed in the introduction to this 
article and in view of the sharply reduced statistics at higher energy and complications from 
the "radiative return" and W pair production we prefer to defer this analysis to a later date. 
Parity violating effects in hadronic collision from SUSY virtual contributions such as those 
treated here in Z decay have also received recent attention [22"]. 
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FIG. 2. Results for different schemes vs universal scalar mass 
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